Strongly Correlated Electrons on a Silicon Surface: Theory of a Mott Insulator 
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We demonstrate theoretically that the electronic ground state of the potassium-covered Si(lll)- 
B surface is a Mott insulator, explicitly contradicting band theory but in good agreement with 
recent experiments. We determine the physical structure by standard density-functional methods, 
and obtain the electronic ground state by exact diagonalization of a many-body Hamiltonian. The 
many-body conductivity reveals a Brinkman-Rice metal-insulator transition at a critical interaction 
strength; the calculated interaction strength is well above this critical value. 
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Transport behavior in crystalline materials is governed 
by the excitation spectrum: insulators have a finite gap 
to excitations while metals have zero-energy excitations. 
Band theory accurately describes this distinction in most 
materials: systems with only filled or empty bands are 
insulating while systems with partially occupied bands 
are metallic. However, the band description may break 
down under circumstances when, roughly speaking, the 
energy cost for forming an extended state exceeds the 
cost for forming a localized state. The resulting ground 
state, which arises from electron-electron interactions 
that band theory cannot describe, is known as a Mott 
insulator (| [| . 

Surfaces provide a potentially fertile environment for 
Mott insulators ||]. Electrons occupying surface states 
may localize more readily than in the bulk, due to two 
significant effects: (1) Atoms at surfaces have lower co- 
ordination than in the bulk, raising the energetic cost for 
electron hopping. (2) Surfaces often undergo reconstruc- 
tions, yielding much larger inter-orbital spacings than in 
the bulk. These effects combine to make surfaces nat- 
ural systems to look for Mott insulating behavior. In 
a recent series of experiments, Weitering et al. |^,^| used 
photoemission and inverse photoemission to demonstrate 
that the K/Si(lll)-(v / 3 x V3)-B surface has a gap at the 
Fermi level. Since this system has an odd number of elec- 
trons per unit cell it must be metallic in a band descrip- 
tion, clearly contradicting the photoemission data. On 
this basis, Weitering et al. hypothesized that this system 
(hereafter K/Si-B) is a Mott insulator. 

In this Letter we explicitly demonstrate, by exact so- 
lution of the appropriate many-body Hamiltonian, that 
the electronic ground state of K/Si-B is indeed a Mott 
insulator. The calculation is in three parts. First, we 
use standard density-functional methods to determine 
the geometrical and electronic structure of this surface 
within the local-density approximation (LDA). Second, 
we map the relevant electronic states onto a many-body 
Hamiltonian, which we then solve on a periodic cluster 
using exact diagonalization techniques. Third, we use the 
resulting many-body ground state to compute the zero- 




FIG. 1. Side and top views of the fully relaxed structure 
of K/Si(lll)-( v / 3 x v^-B. K atoms are gray, Si atoms are 
white, and B atoms are black. The hexagon outlines the unit 
cell, and X's denote Si dangling-bond orbitals. 



frequency conductivity or Drudc weight, D, and then 
show that in the infinite system D — > 0, that is, a metal- 
insulator transition occurs in the thermodynamic limit. 

Boron induces a well-known y3 x v3 reconstruction 
of the clean Si(lll) surface 0. Boron substitutes for 
every third Si atom in the second subsurface layer, and 
the displaced Si assumes an adatom position above the 
boron (see Fig. [[]). The electron in the Si-adatom dan- 
gling bond is transferred subsurface, enabling the B atom 
to participate in four covalent bonds. By this mechanism, 
the surface forms a conventional band insulator, leaving 
each Si adatom with an empty orbital extending away 
from the surface. These orbitals form a triangular lattice 
on the surface. In the experiments of Weitering et al., 
K was then deposited onto this insulating substrate until 
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FIG. 2. LDA band structure of K/Si-B at the K saturation 
coverage of 1/3 ML. The Si-related band crossing the Fermi 
level is half occupied by the K 4s electron. The solid curve 
is the optimized fit of the mean- field Hubbard dispersion, Eq. 
(g), to the LDA eigenvalues. The inset shows the Brillouin 
zone of the fundamental unit cell and of the doubled supercell, 
used to compute t/x/Si-B- 



the saturation coverage was reached. 

To determine the equilibrium structure of K/Si-B, we 
have performed extensive LDA calculations. The calcula- 
tions used a slab geometry with three double layers of Si, 
terminated by H, and a vacuum region equivalent to three 
double layers of Si. Total energies and forces were calcu- 
lated using Hamann and Troullier-Martins pseudopoten- 
tials, and a plane- wave basis with a kinetic-energy cutoff 
of 20 Ry, as implemented in the fhi96md code g. Four 
k-points were used for Brillouin-zone integrations. Full 
structural relaxation was performed on all atoms, except 
those in the bottommost double layer, until the rms force 
was less than 0.05 eV/A. We began by first fully relax- 
ing the surface without K present, and then proceeded 
to determine the equilibrium coverage and geometry of 
the K-saturated surface. 

Experimentally, coverage is monitored via the elec- 
tron work function: at the saturation coverage, the low- 
temperature work function reaches a minimum The 
absolute K coverage is not known from experiment, so 
it must be determined theoretically. We calculated the 
work function for the lowest-energy arrangement of ad- 
sorbates at coverages of 1/6, 1/3, 2/3, and 1 monolayer 
(ML), and find a minimum at 1/3 ML, in agreement 
with the conclusions of Weitering et al. At all cover- 
ages, the experimental photoemission spectra show that 
the Si adatom backbond state persists upon K deposi- 
tion, suggesting the K adsorbates do not break the Si 
adatom bonds |J. We therefore assume that at these 
coverages the adsorbates do not destroy the underlying 
reconstruction. The resulting minimum energy configu- 
ration at the saturation coverage of 1/3 ML is shown in 
Fig. [|. The K adsorbates are in the H3 hollow site, with 
the Si adatom slightly shifted from its position with no 
K present. A metastable state with the K adsorbate in 
the T4 hollow site has an energy 0.1 eV higher per unit 
cell. 



At coverages below 2/3 ML, one expects the K 4s elec- 
trons to partially occupy the surface state arising from 
the empty Si-adatom orbitals. At 1/3 ML there is one 
K per Si orbital, so in the band description the single 
surface band is half occupied; this simple picture is con- 
firmed by the calculated LDA band structure shown in 
Fig. ^. Clearly this system must be metallic within band 
theory. To investigate the importance of electronic inter- 
actions not included in band theory, we derive a single- 
band Hubbard model for the half-filled surface state. The 
Hamiltonian is 



H 



ijcr 



(1) 



where cj CT creates an electron with spin a on site i, and 
nip = ct c icr is the number operator. The sites corre- 
spond to the empty Si orbitals that the K electrons are 
doping. The amplitude for hopping from orbital i to or- 
bital j is given by £y, and £y = tji. There is a Coulomb 
energy cost of U to occupy an orbital with two electrons. 

Our approach to determining the parameters of the 
Hubbard model is similar to other first-principles ap- 
proaches J|,[l(J • We first solve the Hubbard model in the 
mean-field (MF) approximation, and then require that 
the resulting single-particle energies optimally reproduce 
the corresponding LDA spectrum (which is also a mean- 
field theory) throughout the zone. In the MF approxi- 
mation, the up electrons move in the average potential 
generated by the down electrons (and vice- versa) , so the 
MF Hubbard Hamiltonian for the up electrons becomes 



_H" MF = 



(2) 



where (n^) is the average density of the down electrons 
on site i. We assume a paramagnetic state in both the 
LDA and the MF solution to the Hubbard model. 

To determine the hopping amplitudes, Uj, we fit 
the single-particle eigenvalues in the MF Hubbard so- 
lution to the LDA eigenvalues of the surface band at 
100 special fc-points. Note that although the K ad- 
sorbates break the three-fold rotational symmetry of 
the substrate, the LDA band structure remains nearly 
isotropic, and so we assume isotropic hopping. We al- 
low hopping between nearest-neighboring and second- 
nearest-neighboring sites, with amplitudes t\ and t% 
respectively — Third-nearest-neighbor hopping was found 
to be insignificant. Thus the dispersion is given by 

e(k) = 2ti cos(fc y ) -I- 2cos(fc a; -v / 3/2) cos(fc 2/ /2) 

+ 2t 2 cos(k x a/3) + 2 cos(fc a! V / 3/2) cos(3fcj,/2) . (3) 

The optimized amplitudes are t\ = 66 meV and t% = —24 
meV. A plot of the fit and the LDA eigenvalues along 
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FIG. 3. The single-particle energy e+ at the edge of the 
Brillouin zone (the X point), as a function of the electron 
transfer Sn. The squares are LDA supercell results; the solid 
line is the best fit of Eq. (|) to the LDA results for \8n\ < 0.3. 



high symmetry directions is shown in Fig. g; the fit is 
very good, with a rms error of 42 meV. 

To determine the intra-orbital Coulomb repulsion U, 
we subject the system to a density fluctuation by moving 
charge from one Si orbital to another. The optimal inter- 
action parameter, Uk/s'i-b, is then determined by requir- 
ing the LDA solution and the MF Hubbard calculation to 
respond identically. In order to maintain overall charge 
neutrality, a supercell calculation with two Si orbitals is 
required. At the edge of the Brillouin zone, the param- 
agnetic single-particle eigenvalues for a shift of Sn of an 
electron from one orbital to another take the simple form 



e± = ±USn/2, 



(4) 



up to an overall constant which we take to be zero. The 
fit of e + to the LDA eigenvalues is shown in Fig. ||. For 
small charge shifts, the LDA eigenvalues are nearly lin- 
ear, and we obtain Uk/Si-b = 1-23 eV. For larger charge 
shifts, additional bands in the LDA calculation enter that 
are not present in the single-band Hubbard model, and 
the LDA eigenvalues drop very slightly below the MF 
Hubbard eigenvalues. As a check of the reliability of 
this approach, we also determined Uk/si-b by fitting the 
change in the total kinetic energy in the MF Hubbard so- 
lution to the LDA. We find Uk/si-b ~ 1-2 eV, consistent 
with the above result. 

Having determined the parameters of the Hubbard 
model describing K/Si-B, we now solve this model ex- 
actly (using a periodic 16-site cluster) to obtain the 
many-body electronic ground state. The number of 
states in the Hilbert space grows exponentially with the 
number of sites in the cluster: in the Hubbard model, 
each site can have zero, one (either up or down), or two 
electrons, so the number of basis states in the space of an 
iV-site system is 4^. The symmetries of the Hamiltonian 
make the matrix block diagonal, but with 16 sites the size 



of the largest block is still more than 10 7 x 10 7 . Conven- 
tional algorithms obviously cannot diagonalize matrices 
this large, and so we use the Lanczos algorithm to deter- 
mine the exact ground state [ pl| . Storing the Hamilto- 
nian and three basis vectors in memory on an IBM SP2 
required 64 nodes with 1 Gb memory each; the compu- 
tation required about 600 CPU hours. 

To distinguish quantitatively between metallic and in- 
sulating behavior, we calculate the zero-frequency con- 
ductivity or Drude weight, D, of the many-body ground 
state. The Drude weight provides a definitive way to 
distinguish metals from insulators irrespective of the ap- 
plicability of band theory: in the thermodynamic limit, 
metals have non-zero Drude weight, while for insulators 
D = 0. Kohn first showed that D may be calculated from 
the variation of the ground-state energy with respect to 
an applied vector potential [j^,|l2|| , 



D = d 2 E/d(f 2 . 



(5) 



We use this technique to determine the Drude weight of 
our Hubbard cluster as a function of the interaction pa- 
rameter U. The results are shown in Fig. ||. At small U, 
the system is in the metallic regime and so D is large, 
and it decreases monotonically with increasing U . D 
never becomes zero, because only an infinite system can 
undergo a true metal-insulator transition. Instead of a 
critical interaction strength, our 16-site cluster shows a 
transition region in the vicinity of the physically interest- 
ing interaction strengths, near U = 1 eV. 

In general, there may be level crossings in the ground 
state as a function of U. For our hopping amplitudes, 
this is not the case — The ground state evolves adiabati- 
cally with increasing U, and the Drude weight in Fig. |] 
is a continuous function. We see no evidence of the inter- 
mediate (semi-metallic) phase that was found in recent 
slave-boson studies of the Hubbard model on a triangular 
lattice with i 2 =0@. 

To show that K/Si-B lies on the insulating side of a 
Mott-insulator transition, we need the critical interac- 
tion strength for the transition in order to compare with 
our calculated interaction strength, ZTk/Si-b- We can ob- 
tain this value by extending our calculated D(U) from 
the exact finite-cluster result to the thermodynamic in- 
finite limit, using the form derived with the Gutzwiller 
approximation |l|] for the infinite system, 



A>o(t/) oc 



1 - (U/U c f, 

o, 



u <u c 

U>U F . 



(6) 



where U c is the critical interaction strength for the metal- 
insulator transition. This functional form fits our 16-site 
results at small and intermediate values of U extremely 
well, as shown in Fig. From this fit, we estimate the 
critical interaction strength for the Mott metal-insulator 
transition in the infinite system occurs at U c = 0.95 ± 
0.02 eV. Our calculated value for the physical system, 
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FIG. 4. The Drude weight, D, as a function of the inter- 
action parameter U. The circles are exact results from the 
Hubbard model on a 16-site cluster, and the curve is a fit of 
the Hubbard results (for U < 0.7 eV) to the infinite limit 
given by Eq. (§). 



C^K/Si-B = 1-23 eV, is well above this critical value, and 
establishes that K/Si-B is indeed a Mott insulator. 

In the Mott-insulating regime, the localized electrons 
will interact via an antiferromagnetic Heisenberg ex- 
change coupling with strength Jy = 4t?-/Z7 @. The 
Heisenberg model on a triangular lattice is frustrated. 
The nearest-neighbor model has strong three-sublatticc 
correlations, but whether the correlations are long- 
ranged is controversial fig] . In the Hubbard model for 
K/Si-B, the second-neighbor hopping is substantial, so 
the second-neighbor antiferromagnetic Heisenberg cou- 
pling in the Mott-insulating limit will be significant. This 
coupling frustrates the three-sublattice correlations in 
the nearest-neighbor model, so it is unlikely that three- 
sublattice order is established. The ground state will 
likely either establish some other type of collinear order 
fl6| or enter a quantum-disordered regime with no long- 
range order 

To conclude, we have shown that the many-body elec- 
tronic ground state of the K/Si-B surface is a Mott in- 
sulator. Specifically, we have first determined the sur- 
face coverage and morphology of K adsorbed on Si(lll)- 
B using first-principles total-energy methods. We then 
mapped the relevant electronic degrees of freedom onto 
a Hubbard model, which we solved with exact diagonal- 
ization. By calculating the Drude weight of the model, 
we have demonstrated that for the physical parameters 
of K/Si-B, the model has a Mott insulating ground state. 
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